Self— Similar Crossover in Statistical Physics 

V. I. Yukalovf] and S. Gluzman 2 

1 Bogolubov Laboratory of Theoretical Physics 
Joint Institute for Nuclear Research, Dubna 141980, Russia 

2 International Centre of Condensed Matter Physics 
University of Brasilia, Brasilia, DF 70919-970, Brazil 



Abstract 

An analytical method is advanced for constructing interpolation formulae 
for complicated problems of statistical mechanics, in which just a few terms of 
asymptotic expansions are available. The method is based on the self-similar 
approximation theory, being its variant where control functions are defined 
from asymptotic crossover conditions. Several examples from statistical physics 
demonstrate that the suggested method results in rather simple and surprisingly 
accurate formulae. 
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1 Introduction 



In many problems of statistical physics one encounters the so-called crossover phenom- 
ena, when a physical quantity qualitatively changes its behaviour in different domains 
of its variable. To be more precise, we may specify a crossover as follows. Let a func- 
tion f(x) represent a physical quantity of interest, with a variable running through the 
interval X\ < x < x 2 . And let the behaviour of this function, describing some physical 
process, be essentially different near the boundary points x\ and x 2 . Assume that the 
function varies continuously from f(x\) to f(x 2 ), as x changes from x\ to x 2 . Then we 
may say that the function in the interval [xi,x 2 ] undergoes a crossover between f(xi) 
and f(x 2 ). 

Crossover behaviour of different physical quantities is so ubiquitous in nature that 
one could list a plenty of examples. For instance, a number of physical quantities essen- 
tially change their behaviour when passing from the weak-coupling to strong-coupling 
limit [1]. In the theoretical description of crossover there exists a problem which is 
common for practically all physical applications. Real physical systems are usually 
so complicated that describing them equations almost never can be solved exactly. 
However, it is often possible to find asymptotic expansions of solutions in vicinity of 
boundary points. The natural arising problem is how to construct an accurate approx- 
imation for the sought function, valid on the whole domain of its variable, knowing 
only its asymptotic behavior near the boundaries. This problem is aggravated by the 
fact that only a few terms of the asymptotic expansions are usually available. In such 
a case the problem looks unsolvable. 

The most known method of treating the interpolation problem is by using the 
so-called two-point Pade approximants [2-4] or, equivalent to the latter, the Thron 
continued fractions [5,6]. In many cases the two-point Pade approximation yields 
quite reasonable interpolation formulas. However, the usage of this method has not 
become widespread because of the following shortcomings of the Pade approximants: 

(i) When constructing these approximants, one often obtains spurious poles yield- 
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ing unphysical singularities [2-4]. A sequence of Pade approximants may even have 
infinitely many poles [7]. 

(ii) A number of examples are known when Pade approximants are not able to sum 
perturbation series even for small values of an expansion parameter [8,9]. 

(iii) In the majority of cases, except some trivial examples, to reach a reasonable 
accuracy, one needs to have tens of terms in perturbative expansions [4]. While, as is 
emphasized above, in physically interesting problems one often has only a few terms. 

(iv) Defining the two-point Pade approximants one always confronts the ambiguity 
in distributing the coefficients, that is, in deciding which of the coefficients must repro- 
duce the left-side expansion and which the right-side series [1-6]. Such an ambiguity 
increases with the increase of approximants orders, making it difficult to compose two- 
point Pade tables. And for the case of a few terms, this ambiguity makes the two-point 
Pade approximants practically unapplicable. A nice analysis of the latter difficulty was 
done in Ref. [10], where it was shown that, for the same problem, one may construct 
different two-point Pade approximants all having correct left-side and right-side lim- 
its, but differing from each other in the intermediate region by a factor of 40 which 
gives 1000% of uncertainty. This clearly demonstrates that in the case of short series 
the two-point Pade approximation not only cannot provide a reasonable quantitative 
approach but even does not permit to get a qualitative description. The latter concerns 
the general situation, although there can happen some trivial cases when two-point 
approximants make sense even being built with a few perturbative terms. However, 
their application to such few-term cases, in general, is absolutely unreliable. 

(v) The two-point Pade approximants can be used for interpolating between two 
different expansions not always, but only when these two expansions have compatible 
variables [2-4,9]. When these expansions have incompatible variables, the two-point 
Pade approximants cannot be defined in principle. 

(vi) When interpolating between two points, one of which is finite and another is 
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at infinity, one is able to describe at infinity only rational powers [2-4]. The impossi- 
bility to deal with nonrational powers limits the applicability of the two-point Pade 
approximation. 

(vii) The problem of approximating the functions increasing at infinity is especially 
difficult. A two-point Pade approximant can treat only a power-law increase [1-5] 
and is not able to describe other types of behaviour. But in physical problems the 
functions of interest often exhibit at infinity quite different behaviour, for example, 
growing exponentially or following other more complicated ways. In such cases the 
two-point Pade approximants are useless. 

The difficulties listed above are well known and discussed in literature. We have 
cited here only some important references [1-10]. More details on mathematical prob- 
lems in Pade approximation and its applications can be found in several volumes of 
papers, e.g. in Ref. [11]. 

As follows from the above discussion, the two-point Pade approximation in many 
cases is not applicable. It is evident that there is a necessity of developing a more gen- 
eral approach which could overcome the discussed difficulties and would be applicable 
to a larger variety of problems, including those for which the two-point Pade approxi- 
mants cannot be used. It is important that such an approach would provide relatively 
simple analytical formulas for the physical quantities of interest. The advantage of hav- 
ing analytical expressions, as compared to just numbers that could be obtained from a 
numerical procedure, is in the convenience of analysing such expressions with respect 
to various parameters entering into them. Therefore, we keep in mind an analytical, 
rather than numerical, method that would combine relatively simple representations 
for physical quantities with their good accuracy. 

It is worth emphasizing that to derive a new physical formula, valid in the whole 
range of physical variables, is not merely a mathematical challenge but this provides 
new physics, since in the majority of cases realistic physical problems correspond nei- 
ther to weak coupling regime nor to strong coupling limit, but to the intermediate 
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range of parameters. Therefore, it is of great importance for physics to possess a gen- 
eral mathematical tool permitting to derive explicit crossover formulas for arbitrary 
physical phenomena. 

In the present paper we suggest an approach for treating this problem. Our ap- 
proach is based on the self-similar approximation theory [12-22] permitting an accurate 
reconstruction of functions from a few terms of perturbative expansions. The effective- 
ness of the self-similar approximation theory is due to the usage of powerful techniques 
of dynamical theory and optimal control theory. Fast convergence of renormalized se- 
quences is achieved by means of control functions. In the algebraic self-similar renor- 
malization [20-22], we required the algebraic invariance of renormalization-group flow. 
Then, control functions are introduced as powers of a multiplicative algebraic transfor- 
mation. These control functions are defined by the stability and fixed-point conditions 
for a dynamical system called the approximation cascade. In general, the evolution 
equations for a dynamical system can be completed by additional constraints whose 
existence imposes restrictions on the definition of control functions. 

Crossover problem presents an example when additional constraints appear abso- 
lutely naturally. Really, assume that we have a /c-order expansion Pk(x) approximating 
the sought function f(x) in the asymptotic vicinity of the left boundary x — x\. And 
suppose that we are given an asymptotic behavior of this function near the right bound- 
ary x = x 2 - For a moment, take for simplicity that we are given the value f(x 2 ) at the 
right boundary point x = x 2 . When constructing a self-similar approximation f%(x) by 
renormalizing the left boundary expansion Pk(x), we have as an additional constraint 
the right boundary condition fl{x 2 ) = f(x 2 ) . 

We show below that the algebraic self-similar renormalization provides a very con- 
venient tool for treating the crossover problem. This approach permits us to find, 
having just a minimal information about the asymptotic behavior of a function near 
boundary points, a quite accurate approximation for the whole region of the variable. 
In the majority of cases the maximal error of a self-similar approximation is a few per- 
cent and in many cases not more than one percent. In addition to being quite accurate, 
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this approximation is usually given by very simple expressions that are easy to analyze. 
We illustrate the approach by several examples from different branches of statistical 
physics. The variety of considered cases emphasizes the generality of the approach and 
proves that it is a very effective tool for treating arbitrary crossover phenomena. 

Recently, we have applied such an interpolation approach to several quantum- 
mechanical problems [23]. However, what makes the latter principally different from 
the problems of statistical physics is that in quantum mechanics one usually possesses 
quite a number of terms of perturbative expansions, while in statistical physics this 
luxury is rather rare, so that in the majority of cases one is able to derive just a few 
perturbative terms. In the present paper we aim at showing that our interpolation 
method does work for those complicated problems of statistical physics where only a 
few terms of asymptotic expansions are available and other methods are not applica- 
ble. Nevertheless, the self-similar interpolation makes it possible to treat even such 
complicated crossover problems, obtaining simple and accurate formulae. 



2 General approach 

In this section, we formulate the general scheme of the approach not specifying the 
physical nature of a considered function. Let us be interested in a function f(x) of a 
real variable x. Assume that in the vicinity of some point x = x there exist asymptotic 
expansions p k (x, x ), with k — 0, 1, 2, . . ., corresponding to this function, 

f(x) ~p k (x,x ), x -> x . (I) 

Following the algebraic self-similar renormalization procedure [20-22], we define the 
algebraic transform 

P k (x, s, x ) = x s p k (x, x ) , (2) 

where s is yet unknown, and later will play the role of a control function. The transform 
inverse to that in Eq.(2) is 

p k (x, x ) = x~ s P k (x, s, x ) . (3) 
6 



Then we have to construct an approximation cascade with a trajectory bijective to 
the approximation sequence {Pk}- This procedure with all necessary mathematical 
foundations and details has been described in Refs. [13-19]. So, we sketch here only 
the main steps needed for grasping the idea of the method and we concentrate on those 
formulas that permit us to apply the method for crossover phenomena. 
Define an expansion function x = x((p, s, Xo) by the equation 

P Q (x, s,x ) = tp, x = x(tp, s,x Q ) , (4) 

where Po is the first available term from the sequence {Pk}- Introduce a function 

y k (<p, s, x ) = P k {x{if, s, x ), s, x ) . (5) 

The transformation inverse to Eq. (5) reads 

P k (x, s, x ) = y k (P (x, s, x ), s, x ) . (6) 

The family of endomorphisms, {yk}, forms a cascade with the velocity field 

v k ((f, s, x ) = Uk+i(<P, s, x ) - y k (ip, s, x ) . (7) 

The trajectory of the cascade {y k } is, by definitions (5) and (6), bijective to the approx- 
imation sequence {Pk}- Embedding the approximation cascade into an approximation 
flow [16-19] and integrating the corresponding evolution equation, we come to the 
evolution integral 

Jp k v k {f,s,x ) 

in which Pk = Pk{x, s,Xq) is any given term from the approximation sequence 
{Pk}', P k +i = Pk+ii x i s i T i x o) is a self-similar approximation representing a fixed point 
of the approximation cascade; and r is an effective minimal time necessary for reaching 
the fixed point. 

Recall that we started with a sequence {pk} of asymptotic expansions for the con- 
sidered function /(x). Then we passed to the sequence {Pk} by means of the algebraic 
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transformation (2). Now we have to return back employing the inverse transformation 
(3). To this end, we set 



F£(x, s, r, x ) = x s Pk(x, s, r, x ) . 



(9) 



The quantities s and r are the control functions guarantying the stability of the 
method, that is, the convergence of the procedure. These functions are to be defined 
by the stability conditions, such as the minimum of multiplier moduli, together with 
additional constraints, like, e.g., boundary conditions. Let us find from such condi- 
tions s = Sk and r = t\~. Substituting these into Eq.(9), we obtain the self-similar 
approximation 



for the function f(x). We retain here the notation for the point xq in order to show that 
the approximation (10) has been obtained by renormalizing p k (x, xo) which, according 
to Eq. (1), is an asymptotic expansion of f(x) in the vicinity of the point x = x . 

Now assume that the variable x changes in the interval X\ < x < x 2 and that the 
asymptotic behavior of a function f(x) is known near the boundaries of this interval. 
The latter means that in Eq. (1) we have to put, instead of x , either x\ or x 2 . Let 
us take, for concreteness, the boundary points x\ — and x 2 — > oo. Then we have 
two types of expansions, pk(x,0) and p k (x,oo). Following the procedure described 
above, we can construct, in the place of Eq. (9), two quantities, F£(x,s,t,0) and 
Fk(x,s,r, oo). 

As is discussed above, the control functions s and r are to be defined from stability 
conditions plus additional constraints. The natural such constraints for the crossover 
problem can be formulated as follows. Suppose we have constructed the renormalized 
expression F£(x,s,t,0) starting from the left asymptotic expansion p k (x,0). By this 
construction, the function F£ has correct asymptotic behavior near the left boundary. 
But in order to correctly represent the sought function in the whole interval of x G 
[0, oo), the renormalized expression F£ must have the correct asymptotic behavior when 



ft(x, Xq) = F£(X, S k , T fc , X ) 



(10) 
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approaching the right limit. This implies the validity of the condition 

Jim \F*(x, s, r, 0) - Pi (x, oo)| = , (11) 

imposing constraints on s = s k and r = T k . We shall call Eq. (11) providing the correct 
crossover behavior from the left to the right boundary the left crossover condition. The 
quantity s k can be called the left crossover index, and r k , the left crossover time. For 
the self-similar approximation (10) we get, in this way, 

r k (x,0) = F*(x,s k ,r k ,0), (12) 

which may be named the left self-similar approximation, or the left crossover approx- 
imation. 

The analogous reasoning works, as is clear, when we are considering the crossover 
from the right to left. Then we obtain the right crossover condition 

\im\F k (x,s,r,oo) -pj(x,0)\ = , (13) 

imposing constraints on s = s k and r = r k , thus defining the right crossover index 
s k and the right crossover time r k . As a result, we come to the right self-similar 
approximation, or the right crossover approximation 

f k (x, oo) = F k (x, s k , r k , oo) . (14) 

In general, from the left and the right approximations, (12) and, respectively, (14), we 
can compose the average self-similar approximation, or average crossover approxima- 
tion 

r k {x) = l -\r k {x,Q) + r k {x,^)] . (is) 

The suggested general approach to reconstructing crossover functions can be em- 
ployed for any crossover phenomena. In particular applications, it can happen that we 
possess a reliable asymptotic expansion only from one side of the crossover domain, 
and from another side just one term is available. In this case, as is clear, we are not 
able to construct both left and right self-similar approximations, but only one of them. 
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Nevertheless, such one-side approximations are usually quite accurate, as we show by 
examples in the following sections. The possibility of constructing accurate approxi- 
mations, when we have a perturbative series only from one side of the crossover region 
and a sole asymptotic term from another side, is very important since this situation 
constantly occurs in realistic physical problems. We shall demonstrate in what follows 
how it is possible to improve the accuracy of such one-side approximations by com- 
bining the terms of a given one-side series and defining the crossover indices so that 
to satisfy the asymptotic limit from another side, in accordance with the crossover 
conditions (11) or (13). 

In order to emphasize that the suggested approach does work even for the cases 
with a very scarce information about the sought function, let us consider a simple 
example. Suppose that we know the asymptotic behavior of a function near the left 
boundary, where x — > 0, only in the linear approximation 



as x — > oo. In such an extreme case of minimal information, it looks like there is no 
a regular way of recovering the function for the whole axis < x < oo. However, 
our approach, based on the idea of self-similarity, permits us to recover the sought 
function. 

Following the procedure described above, in the place of Eq.(9), starting from ex- 
pansion (16), we obtain 





(16) 



And assume that only one asymptotic term is known from the right side, 



Pi(x, oo) 



A, n ^ , 



(17) 




(18) 



With Eqs.(17) and (18), the left crossover condition (11) reads 
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From here the left crossover index si and the left crossover time T\ are uniquely defined 
as follows: 



a 



si = -n , n= n— — J . (19) 

Substituting these values into Eq. (18), as is prescribed by Eq. (12), we recover the 
left crossover approximation 



/i*(x,0) = a 



1 f ( — x 

Oq J 



-i n 



(20) 



At large rr — > oo, expression (20) reduces to the limit (17). When x — > 0, then we have 
the linear behavior 

fi(x,0) ~ a + ajx , 

where = na (A/a ) 1//n is the renormalized coefficient. Such a renormalization is 
typical of renormalization group techniques, as is discussed in Refs. [20-22]. 

Thus, even having so scanty information about the asymptotic properties of a func- 
tion, as in the above example, our approach allows us to reconstruct, in a systematic 
way, the function for the whole domain of its variable. This reconstruction becomes 
possible owing to the idea of self-similarity which our approach is based on and due to 
the convenient introduction of control functions through the algebraic transformation. 
The idea of self-similarity, complimented by the property of algebraic invariance, elimi- 
nates the umbiguity typical of divergent series in the standard perturbative approaches. 
In the sections that follow, it will be shown that the accuracy of so constructed self- 
similar approximations is rather good. 

Note that achieving good accuracy with a limited number of terms of an asymptotic 
expansion should not be treated as surprising. Asymptotic series are known to provide 
reasonable accuracy when up to some optimal number of terms are taken [24], the 
subsequent terms only spoil the picture being, in this sense, excessive. Whether there 
are such excessive terms or not is decided, in our approach, by stability and crossover 
conditions. As soon as these are satisfied, there are no excessive terms. And if adding 
more terms does not allow us to satisfy these conditions, the added terms are to be 
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considered excessive. Fortunately, the real life and realistic physical problems are so 
complicated that we practically never have excessive terms, but vice versa, have to deal 
with very short expansions containing only a few terms. 



3 Zero— dimensional model 

For illustrative purpose, we start with a simple model example. Consider the partition 
function of a zero-dimensional anharmonic model, represented by the integral 

J(g) = —= / exp (-X 2 - gx 4 ) dx, (21) 



with the integrand possessing a single "vacuum" state, located at the point x — 0. 
The weak-coupling expansion of this integral in powers of the coupling parameter g, 
around the vacuum state, leads to divergent series, 

J( g )~a + bg + cg 2 + ... , (g - 0), (22) 

where 

i h 3 105 

a—i, o — — - , c = — — . 

4 ' 32 

The so-called strong-coupling expansion, in inverse powers of the coupling constant, 
can be written down as well: 

J(g) ~ Ag~ 1/4 + Bg^ + Cg^ 4 + . . . , (g - oo) , (23) 

with 

_ 1.813 _ 0.612 _ 0.227 

V 71 " v 71 " v 71 " 

Following the approach of Section 2, one can derive the right crossover approximation, 

r(g,oo) = aA(A 2 + aV^y 1/ \ (24) 

with the right crossover index s = 1/2 and crossover time r = — A 3 /(2a 2 B) = 1.55. 
At # = 1, the percentage error of formula (24), is equal to —7.38%, while the maximal 
error is reached at g — 0.35 and equals —7.96%. 
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The left crossover approximation is given as follows: 

J*(g,0) = aA(A A + a i g y 1/ \ (25) 

with the left crossover index s = 1/4, and crossover time r = — a 5 /(4A A b) = 0.304. At 
g — 1, the percentage error of Eq. (25) is 10.13%, while the maximal error at g — 2.5, 
is equal to 10.53%. We conclude, that the crossover approximations (24) and (25) may 
be viewed, correspondingly, as the lower and upper bounds for the integral (21). The 
average, defined by Eq. (15), 

J*(g,0) + J*(g, oo) 

J (g) = 2 ' 

possesses the correct leading asymptotes and approximates the exact result at g — 1 
with the percentage error of 1.37%. And the maximal error, at g — 3, is 2.21%. 



4 Lattice gauge model 

The vacuum energy density f of the (3+l)-dimensional SU(2) lattice gauge model 
in its weak-coupling, asymptotically free regime, may be presented in the form of an 
expansion in powers of the parameter x = 4/g 4 , where g stands for the coupling [25]: 

fo ~ Ax + + ... , (x -> oo) , (26) 

where 

A = -6, B — 7.1628 . 
In its strong-coupling limit, fo can be presented as follows [25]: 

fo ~ ax 2 + bx 4 + . . . , (x -> 0) , (27) 

with 

a = -1 , b = 0.03525 . 

Because of the interfering roughening transition, the quality of the high-order terms in 
the strong-coupling expansion is doubtful [25], so we use only its leading terms. The 
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left crossover approximation can be readily written down as 

-1/2 



where we have used 



f*(x, 0) = ax 2 



1 

2 



1 + 



a V 

a) 



x 



(28) 



2bA 2 



The numbers generated by formula (28), practically coincide in the region x e [0.1, 1.1] 
with estimates obtained in [25] from the strong-coupling approximants. The right 
crossover approximation can be obtained as well, but its accuracy is worse than that 
ofEq. (28). 



5 One-dimensional Bose system 

The ground-state energy of the one-dimensional Bose system with the (^-functional 
repulsive interaction potential is known in a numerical form from the Lieb-Liniger 
exact solution [26]. It is desirable, nevertheless, to have a compact analytical expression 
for the ground-state energy e(g) as a function of the 5-function strength g, valid for 
arbitrary g. In the weak-coupling limit an exact analytical result is known: 

e(g) ~g, (g - 0) , (29) 

obtained in Ref. [26], while in the strong-coupling limit another exact result, obtained 
by Girardeau [27], is available: 

e(g) ^ y , (g - oo) . (30) 

The higher-order terms in these expansions were derived by approximate methods, the 
next term in the weak-coupling limit being ps bg 3 ^ 2 , and in the strong-coupling limit 
Bg~ l . We shall not use the approximate values for the coefficients, b and B (see e.g. 
[13] and references therein), writing instead trial expansions and determining the coef- 
ficients by matching the two asymptotic forms for the ground state energy. Following 
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the standard approach of Section 2, we obtain the right crossover approximation: 



with the right crossover index s = 1 and B = (ir 2 /3) 2 . Although, Eq. (31) can 
be further simplified, we leave it in present form in order to stress the origin of its 
different parts. Simple expression (31) works with surprising accuracy of about 1 — 2%, 
up to g ~ 10, till there are numerical data available for comparison [26]. The left 
crossover-type expression can be written as well, following the standard procedure, 
but its accuracy is inferior to that of Eq. (31). 

6 One— dimensional ferromagnet 

Low-dimensional magnetic systems give a plenty of examples of the crossover phe- 
nomena, when only the asymptotic behavior with respect to different parameters, such 
as spin, temperature etc, is known and the intermediate region, in most of the cases, 
could be reached only numerically. The crossover self-similar approximations offer 
simple analytical expressions for the intermediate region. We put below, for simplicity, 
the exchange integral J = 1. 

6.1 Zero— field thermodynamics 

The free energy F and magnetic susceptibility x °f the one-dimensional Heisenberg 
ferromagnet of spin S, within the framework of the spin-wave approximation, valid at 
temperatures T — > 0, has the form of an expansion in powers of T [28]: 




(31) 



F ~ a(S)T 3 / 2 + b(S)T 2 + . . . , (T — > 0) , 



(32) 



in which 





1 



and 



X ^ A(S)T~ 2 + B(S)T 3/2 , (T 



0), 



(33) 
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where 



A(S) = ^S 4 , B(S) = -A(S) 



3C(l/2) 



3 ' vw v ~' V2ttS 

As T — > oo, a different asymptotic behavior happens [29]: 



F ~ -Tln(l + 25) (T^oo) 



and 



X ^ 



4g(g+l) 
3T 



(T -> oo). 



(34) 



(35) 



Applying the standard approach of Section 2, we obtain for the free energy the following 
left crossover approximation, corresponding to the left crossover index s — 1, 

T 3 / 2 a 2 (S 1 ) 



F*(T, 5) = a(S): 



T = 



1 - [b(S)/a(S)]rT 1 / 2 ' ' ~ b(S) ln(l + 25) 
and the expression for specific heat C* = —Td 2 F*(T, S)/dT 2 as 

r*(r^- ±m n 3 ( o, -MS) + b(S)TT^ 
6 (T ' S) ~ ~l T a {S) l-a(S) + 6(S)tTV2]3 • 



(36) 



(37) 



The position, height and spin-dependence of the maximum occurring in the expression 
for C*(T, 5) are in qualitative agreement with numerical results for finite chains [30]. 
The left crossover approximation for the renormalized susceptibility is 



X 



MS) 

J^2 



1 + 



B (s)_ rTl/2 



2A(S) 



MS) 
B(S) 



165(1 + 5) 
3A(S) 



1/2 



(38) 



with the left crossover index s = — 2. The expressions (36) and (38) are very accurate 
for 5 = 1/2, where they practically coincide with the results of a numerical solution of 
the thermodynamic Bethe-ansatz equations [28]. 



6.2 Spin waves at finite temperatures 

Variational theory, as applied at low temperatures [31], gives the temperature- 
dependent expression for the spin-wave energy Uk for 5 = 1/2 in the form 



uj k = 2Z(T)\sm(k)\ , Z(T)~-n 



j 2 /T x 2 
~ 3 \2 



(T->0), 



(39) 
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being at T = completely in agreement with the exact results [28]. In order to find 
the behavior of Z(T) at arbitrary T, we continue it from the region of T — > self- 
similarly, along the most stable trajectory, with the crossover index s, determined by 
the condition of the minimum of the multiplier [20-22] 



m(T, s) = l--T 



from where 



-i 



s{T) = -T 2 (l- 1 ^) T<VQ, 



s -> oo , T>VQ. 
This gives the left crossover approximation 

^^(mtWfT^ ' T<V ~ 6 ' (40) 

Z*(T) = i7rexp(-T 2 /6) , T > V6 . (41) 

Formulae (40) and (41) suggest that the spin waves should survive at least up to 
T ~ \/6, and become exponentially "soft" above this temperature. Note that in this 
particular case the left self-similar approximation plausibly reconstructed the function 
for arbitrary temperatures, even not knowing beforehand the asymptotic behavior at 
T — > oo. 



6.3 Field— dependent part of free energy 

It is believed that the magnetic field-dependent part of the free energy of the one- 
dimensional Heisenberg ferromagnet is independent on spin and scales as p = h/T 2 
(h denotes the magnetic field), with the scaling function independent on the value of 
spin [32]. For the classical ferromagnet, both low and high field behavior of the field 
dependent part of the free energy 5F(p) are known [32-35] in the simple form: 

T- 2 5F(p)~ap 2 + bp\ a = ~, b=^ (p « 1), (42) 
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T~ 2 5F{p) ~ Ap + Bp 1 / 2 , A = -l, B = l (p»l). 



(43) 



The left crossover approximation is controlled by the crossover index s = 1/2 and 
crossover time r = — a 3 /(2A 2 B) = 0.227, yielding: 



T- 2 5F*{p,0) = ap 2 
while the right crossover approximation is given by 

T' 2 5F*(p, oo) = Ap 2 

with 



-1/2 



(44) 



, A\ 1 / 2 

p"M^) 



-2 



(45) 



A / Ax 1 / 2 

—40 = iM - 



s = 2, 

Both expressions (44) and (45) are in good agreement with the known results [32-35]. 



7 Flexible polymer coil 

The calculation of the so-called expansion function a 2 (z) of a flexible polymer coil is 
of long standing interest in polymer science [36-38]. This quantity defines the ratio of 
the mean square end-to-end distance < R 2 > of the chain to its unperturbed value 
< R 2 > = Nl 2 , where N is the number of segments with the length / each, so that Nl 
is the contour length of the chain, 

< R 2 > 



a{z) ~<R 2 > ' 
as a function of a dimensionless interaction parameter z. The latter is 

BN 



z = 



7t/ 2 



(D = 2) 



(46) 



(47) 



for the two-dimensional case and 



: = I — I —^r- (D = 3) 



,2,7 (3 - -, (48) 

for the three-dimensional coil, where B is the effective binary cluster integral for a pair 
of segments. 
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When the excluded volume interaction is very weak, a perturbation theory leads 
[39] to an asymptotic series 

a 2 {z) ~1 + Y.anz n (*->0), (49) 

71=1 

in which the coefficients for the two-dimensional case are 

ai = - , a 2 = -0.121545 , a 3 = 0.026631 , a 4 = -0.132236 (D = 2) , 

and for the three-dimensional coil they take the values 
4 

ai = - , a 2 = -2.075385 , a 3 = 6.296880 , a 4 = -25.057251 , 

a 5 = 116.134785 , a 6 = -594.71663 (D = 3) . 
The asymptotic result for the strong coupling limit [40] is 

a 2 ^) ~ A x z p + A 2 z 7 (z -> oo) . (50) 

Using our method of self-similar interpolation, we obtain from (49) and (50) 

o£(z) = (l + A\" i zy (51) 

in the first approximation. The second approximation gives 

[1 + C lZ ) 2 -^ + C 2 z 2 ] m , (52) 



oit(z) = 



where 



2A 2 \ 1/(2_/3+7) 

^1 = I " ,1-27/3 ) ' C 2 = -Ax . 



These formulae can serve for both the two- as well as for three-dimensional coils. We 
shall concentrate on the latter case for which accurate numerical data for a 2 (z) are 
available [40] in the whole range of z G [0, oo). Then in the strong coupling limit (50), 
one has 

A x = 1.5310 , A 2 = 0.1843 , (3 = 0.3544 , 7 = -0.5756 (D = 3) . 

(53) 
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The coefficients in (52) become 



d = 6.5866 , 



C 2 = 11.0631 , 



2-/3 + 7 = 1.07. 



In this way, from (52) we obtain 



al(z) = (1 + 6.5866^) L07 + 11.0631r 



2 1 0.1772 



(54) 



The self-similar approximation (54) is accurate, within 0.4% of error, in the full range 
z > 0, as compared to numerical calculations [40]. This formula (54) practically 
coincides with the phenomenological extrapolation expression 



obtained by Muthukumar and Nickel [40] by means of a fit to numerical data. 

In conclusion, we have developed the method of self-similar interpolation for de- 
riving explicit interpolation formulae for difficult crossover problems of statistical me- 
chanics. This method, as is illustrated by several examples, is general, simple, and 
accurate. 
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